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We present a program for the evaluation of full unpolarized cross sections for the W ± and Z° production 
in the narrow width approximation at NNLO in perturbative QCD using Mellin space techniques. 
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1 Introduction 

The Drell-Yan process, originally described in the context of the parton model [ 1 1, concerns the production 
of a lepton pair of large invariant mass in hadron-hadron collisions. With the increase of the centre of mass 
energy at particle accelerators, the Drell-Yan process led to the discovery of W and Z a bosons at UA1 
and UA2 experiments J2]|3]. Since then the properties of massive vector bosons have been studied in great 
detail. At present the production of W and Z° provides an important benchmark for the LHC and a test 
of the Standard Model (SM) in a new range of centre of mass energies J4). 

As guaranteed by the factorisation theorem Q, one can separate the physics of soft energy scales from 
the physics at hard energy scales where perturbation theory applies. The higher order QCD corrections 
to the Drell-Yan process have been calculated up to next-to-next-to-leading order (NNLO), see [6-8 1 and 
references therein. The full cross section is obtained as a convolution with the parton distribution functions 
(PDFs) that encode the non-perturbative information. 

In this paper, we present a program for evaluation of the full inclusive cross section for W and Z° 
production in a fast and accurate way using a Mellin space approach. After a brief description of the basic 
ingredients of the calculation we give formulae for the Mellin transforms. We then present a comparison 
with the code ZWPROD (TJtEl and discuss possible applications and extensions within this framework. 



2 Formalism 

We consider the inclusive production of a single vector boson V = W + , W~ or Z° in hadron-hadron 
collision with a centre of mass energy s which subsequently decays into a lepton pair of an invariant mass 
Q 2 . The decay of the vector boson is treated within the narrow width approximation which replaces the 

this process can be expressed as 



propagator by a delta function such that Q 2 = M v . We consider massless quarks. The cross section for 



a tlltl ^ v ^ hl2 (s)=xa v ^ llh W v (x,Q 2 ), x = Q 2 /s, (1) 

where (j v ^ 1112 represents the kinematically independent part of the Bom level subprocess qq — > V — > lih 
(the point-like cross section) multiplied by the appropriate branching ratio. The exact form of the point 
like cross section can be found in Ref. 0, formulae (A. 10) and (A. 1 1). Q 



* Corresponding author E-mail: petra.kovacikova@desy.de, Phone: +00493 376277 193 
' There is an extra factor of 2 in the denominator of the formula (A. 11), corrected by 1 9 1 
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The structure function W (x, Q 2 ) is written as a convolution of two parton distribution functions 
f a (^ 2 f) and fb(fif) and a hard scattering cross section represented by coefficient functions A ab , 

W V (X,Q 2 ) = C Zb[fa(t4)®fb($)®Kb(Q 2 ,»%V 2 r )](x). (2) 

a,b=q,q,g 

The perturbative coefficients are known up to NNLO []7]|8], 

A^°(x, Q 2 , »}, ^ = E ^Af) (*, Q 2 , „» , £). (3) 

The factor C a v b in Eq.© contains information about couplings of vector bosons to partons a and b. For 
the detailed form of the Eq.© we refer the reader to the paper of Hamberg, Matsuura and van Neerven J7] 
whose notation we follow closely 0. The convolution sign represents an integral 

Cfi ®/ a ® ■•■®/fe) 0) 

dx-i / dx 2 ... dx k 5(x - X!X 2 ■ ■ ■Xk)fi(x 1 )f 2 (x 2 ) ■ ■ ■ f{xk)- (4) 



In principle one can perform the integrals in Eq.(|2} directly however, the problem is much better addressed 
after transforming to Mellin space, 

/(TV) = f 1 dxx N -\f(x), (5) 
Jo 

This transformation turns the integrals in Eq.(|2) into ordinary products such that the structure function 
reads 

W V (N,Q 2 ) = J2 Cl b f a (N,Q 2 )f b (N,Q 2 )A ab (N,Q 2 ) M/ = /V = Q 2 , (6) 

a,b=q,q,g 

and therefore it is possible to evaluate it in a fast and efficient way. The formula for the inverse Mellin 
transform defines how to recover the original momentum space result, 

f c+ioo 

(7) 



W v (x,Q 2 ) = 7— dNx- N W v (N,Q 2 ), 

ZTTl Jc—ioo 



where c represents a point on the real axis such that all poles TV, in the function W(N, Q 2 ) lie to the left 
from c. Further on, we will refer to functions in Mellin space as N space functions and functions from 
momentum space as x space functions. 



3 Implementation 

The main ingredients of the calculation are the coefficient functions up to NNLO and the parton distribution 
functions in Mellin space in terms of a complex variable N. The condition N e C is required for the 
numerical evaluation of the inversion formula (|7}. For this we adopted the technique implemented in 
QCD-PEGASUS iflOl . The complex integral (0 is rewritten in terms of an integral over a real variable z 

1 f°° 

W(x,Q 2 ) = - dzlm[e ltp x- c - ze ^W(N,Q 2 )} JV = c + zexp^eC (8) 
* Jo 

and evaluated using Gaussian quadratures. The parameter <f> > ir/2 represents the angle with respect to the 
positive real axis. Since the rightmost pole of the structure function is N max = 1, we chose c = 1.5. These 

^ Several typos appearing in Ref. (Jj have been pointed out in (9) 
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values as well as the maximum value of the integration variable z are flexible and can be modified by the 
user in the main program if desired. For a more detailed description of the shape of the integration contour 
we refer to the QCD-PEGASUS manual ifTOl . 

The coefficient functions in N space were published previously in Ref. iTTTI . including corrections to 
the previous literature. The corresponding FORTRAN code is DY . f used together with ANCONT fPH . We 
performed the Mellin transforms starting from the x space expressions lfl"3l using the harmpol pack- 
age fl4l . The results can be expressed mostly in terms of complex-valued simple harmonic sums 11511161 
and several more complicated ones which we approximated by using the minimax methoc0 worked out in 
detail in lfl2ll previously^ The absolute accuracy of our approximation is better than 10~ 9 over the whole 
kinematic range. 

At the moment there are two options for the input parton distribution functions in N space. A toy 
input corresponds to the one used for the 2001/2 benchmark tables [1811 and is used for comparisons with 
ZWPROD (TJISl assuming no evolution of PDFs. The general form reads 

xfi,toy( x > Po) = nx a {l - x) b 7 i = q,q,9, n,a,beR, (9) 
which is in Mellin space represented by an Euler beta function 

f ittoy (N,iJ, 2 )=nj3{a + N,b + l). (10) 

The second option for the PDF input is using the FORTRAN code QCD-PEGASUS QUI which can be 
linked to our program. 
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Fig. 1 Cross section for W~ production up to NNLO in the narrow width approximation using the toy 
parton distribution functions and a fixed value of the strong coupling constant. Upper part: The full cross 
section. Lower part: Relative accuracy with respect to the ZWPROD. 



4 Results And Outlook 

There are several programs on the market using the standard momentum space evaluation [19-21] which 
can provide a cross-check for our N space calculation. We performed comparisons of the full cross sec- 
tions with a program ZWPROD written by the authors of the original calculation of the NNLO Drell-Yan 



We used the MINIMAX routine implemented in Maple 
Exact expressions were given in 1171 . 
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coefficient functions (7JIS1 ■ The Fig. ([TJ shows a comparison for the W~ cross section using toy input for 
PDFs corresponding to the Eq.rtTOt with no evolution and a fixed value of the coupling constant a s = 0.1. 
The relative accuracy is better than 6 x 10~ 6 in the relevant kinematical range x £ (10~ 4 , 0.8). As an 
intermediate check, we compared the Mellin inversion of N space coefficient functions against the x space 
expressions using a program of Gehrmann and Remiddi EZl for the numerical evaluation of harmonic 
polylogarithms. The framework presented here is suitable for a further implementation of those cross 
sections where A~ space coefficient functions are also available, like Higgs production and deep inelastic 
scattering (DIS) |TTll23]-[26l. The setup is well suited for merging the program with threshold resum- 
mation calculations which are typically performed in Mellin space (see e.g. 11271 ). For the extraction of 
PDFs from and Z production it would be desirable to have an access to the rapidity distributions in 
which case one will need to apply double Mellin transforms of two variables N\ and Af 2 however, this is 
a subject to further study. On the side of PDFs we aim for a direct interface to the LHAPDF grids l28l 
which will allow the user to freely choose any particular PDF set provided within this framework. Recent 
results 11291 on A~ space input parametrizations also allow for more flexible input PDF parametrisations in 
QCD-PEGASUS. Further improvements with respect to the speed of the code are foreseen and together 
with an upgrade on the input PDFs this code can become a tool for PDF fits, where fast and accurate 
evaluations of cross sections are needed. The current version of the C++ code can be downloaded from 
|http : / / www- zeuthen . desy . de/ ~kpetra/ sbp| 
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